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Abstract 

An exact analytic solution has been obtained for a uniformly expanding, neutral, infinitely conducting plasma cylinder in an external 
uniform and constant magnetic field. The electrodynamical aspects related to the emission and transformation of energy have been 
considered as well. The results obtained can be used in analysing the recent experimental and simulation data. 
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Many processes in physics involve boundary surfaces which 
requires the solution of boundary and initial value problems. 
The introduction of a moving boundary into the physics usually 
precludes the achievement of an exact analytic solution of the 
problem and recourse to approximation methods is required UJ- 
yD (see also references therein). In the case of a moving plane 
boundary a time-dependent translation of the embedding space 
immobilizes the boundary at the expense of the increased com- 
plexity of the differential equation. It is the aim of this work 
to present an example of a soluble moving boundary and initial 
Value problem in cyhndrical geometry. 

. The problems with moving boundary arise in many area of 
physics. One important example is sudden expansion of hot 
plasma with a sharp boundary in an external magnetic field 
which is particularly of interest for many astrophysical and lab- 
oratory applications (see, e.g., |4] and references therein). Such 
kind of processes arise during the dynamics of solar flares and 
flow of the solar wind around the earth's magnetosphere, in ac- 
tive experiments with plasma clouds in space, and in the course 
of interpreting a number of astrophysical observations Js Mllll . 
' To study the radial dynamics and evolution of the initially 
spherical or cylindrical plasma cloud both analytical and nu- 
merical approaches were developed (see, e.g., Refs. 11 31 41411 and 
references therein). The plasma cloud is shielded from the pen- 
etration of the external magnetic field by means of surface cur- 
rents circulating inside the thin layer on the plasma boundary. 
Ponderomotive forces resulting from interaction of these cur- 
rents with the magnetic field would act on the plasma surface 
as if there were magnetic pressure applied from outside. After 
some period of accelerated motion, plasma gets decelerated as 
a result of this external magnetic pressure acting inward. The 
plasma has been considered as a highly conducting media with 
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zero magnetic field inside. From the point of view of electro- 
dynamics it is similar to the expansion of a superconductor in a 
magnetic field. An exact analytic solution for a uniformly ex- 
panding, highly conducting plasma sphere in an external uni- 
form and constant magnetic field has been obtained in il2ll . 
The non-relativistic limit of this theory has been used by Raizer 
il3ll to analyse the energy balance (energy emission and trans- 
formation) during the plasma expansion. The similar problem 
has been considered in Ref. |8] within one-dimensional geom- 
etry for a plasma layer In our previous paper ll4ll we obtained 
an exact analytic solution for the uniform relativistic expansion 
of the highly conducting plasma sphere in the presence of a 
dipole magnetic field. In the present paper we study the uni- 
form expansion of the highly conducting plasma cylinder in the 
presence of a constant magnetic field. For this geometry we 
found again an exact analytical solution which can be used in 
analysin g th e recent experimental and simulation data (see, e.g., 
Refs. ll4l llin and references therein). 

2. Moving boundary problem 



We consider the moving boundary problem of the highly con- 
ducting plasma cylinder expansion in the vacuum. Consider a 
cylindrical region of space with radius p - R{t) at the time t 
containing a neutral infinitely conducting plasma which has ex- 
panded at f = (with R{Q) - 0) to its present state from a linear 
source located atp = 0. We assume that at any time t the plasma 
cylinder is unbounded in z direction (i.e. the cylinder is located 
at -oo < z < co). To solve the boundary problem we introduce 
the cylindrical coordinate system (p, (/?, z) with the z-axis along 
the plasma cylinder symmetry axis and the azimuthal angle (f 
is counted from the plane (jcz-plane) containing the vector of 
the constant and homogeneous magnetic field Hq. The angle 
§ between the vector Ho and the z-axis is arbitrary. The un- 
perturbed magnetic field is expressed by the vector potential. 
Ho = V X Ao, where this potential is Ao - 4 [Ho x r] with the 
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cylindrical components 



z p 

A(v = - (yo^oii - Z^o± cos ip) . 



(1) 



The components of the vector Hq are //op = ^oj. cos ip, Hq^ - 
-Hqx simp, and Hq^ = Hq\\. Here i/o± = //osinj? and i/o|| = 
//() cos )? are the components of the magnetic field transverse 
and parallel to the z-axis, respectively. 

As the cylindrical plasma cloud expands it both perturbs the 
external magnetic field and generates an electric field. Within 
the cylindrical plasma region there is neither an electric field 
nor a magnetic field. We shall obtain an analytic solution of the 
electromagnetic field configuration. 

We consider the case of the uniform expansion of the plasma 
cylinder R(t) - vt with a constant expansion velocity v. This 
special case of the uniform expansion falls within the conical 
flow techniques which has been applied previously in Ref. 1 12]. 
From symmetry considerations one seeks a solution for the total 
(i.e., the unperturbed potential Aq plus the induced one) vector 
potential of the form A(r, t) ~ p''a(f ) with ^ = p/cf, where c is 
the velocity of light. Since the external region of the conducting 
cylinder is devoid of free charge density, a suitable gauge allows 
the electric and magnetic fields to be derived from the vector 
potential A. Having in mind the symmetry of the unperturbed 
magnetic field and the fact that the electromagnetic fields do not 
depend on the coordinate z it is sufficient to choose the vector 
potential in the form Ap (r, t) - 0, 



i,(r,0 = ^p;r(f), 



A^(r,f) = //o±p;r(^)sin^, (2) 



where xiO is some unknown function. The components of the 
electromagnetic field are expressed by this function as 



-H^^smip[(x{0]' , (3) 



Hp = Ho^x iO cos (^, H^ 
H,^H^,^\^x'iO+xiO 
E^ = ^^''x' (0 , E, = Ho^i^x' (0 sin 'P, 



(4) 



Ep = 0, where the prime indicates the derivative with respect 
to the argument. The equation for the vector potential A(r, t) is 
obtained from the Maxwell's equations which for the unknown 
function ;\f(^) yields an ordinary differential equation 



x"(0+ ,], \^ x'(0-Q- 



(5) 



This equation is to be solved in the external region p > R{t) 
subject to the boundary and initial conditions. Here R{t) is the 
plasma cylinder radius at the time t. The initial conditions are 
atf = 

A(r,0).Ao, ^^.0. (6) 

at 

The first initial condition states that the initial value of A(r, t) 
is that of a homogeneous unperturbed magnetic field. The sec- 
ond initial condition states that there is no initial electric field. 



Boundary conditions should be imposed at the cylindrical sur- 
face p - R(t) and at infinity. Because of the finite propaga- 
tion velocity of the perturbed electromagnetic field the mag- 
netic field at infinity will remain undisturbed for all finite times. 
Further, no incoming wave-type solutions are permitted. Thus, 
for all finite times H(r, f) ^ Hq at p ^ ca. The boundary con- 
dition at the expanding cylindrical surface is Hp = which is 
equivalent to the relation ;^'(/3) = (see Eq. (O) withyS - v/c. In 
addition imposing that H(r, t) = Ho at p ^ cf we obtain another 
boundary condition xW = 1 ■ 

The solution of Eq. ^ subject to the initial and boundary 
conditions may be written as 
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The complete solution may finally be written in the form at 
vf ^p^cf(orySs;^s^ 1) 

no 
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H(r,0 = Ho, E(r,f) = at p > cf (or ^ > 1) and H(r,0 = 
E(r, = at p < R(f) = vt (or ( < /?). From Eqs. ©, ©, Q 
and ([Tol l it can be easily checked that the boundary condition on 
tiie moving surface, E(/?) = -^[v x n(R)] (or E^{R) = pH^{R), 
E^{R) - -/3H^{R)), is satisfied automatically. 

It is also imperative to determine the surface current den- 
sity induced due to the moving boundary. This can be done 
employing the Maxwell's equation j = ^ (cV x H - ^|. It 
is clear that the surface current has only two components and 
jo - iaS{p - R) with a - z,ip, where ia is the linear surface cur- 
rent. Using Eqs. (|9]l and ( fTOl l as well as the Maxwell's equation 
we obtain 



i, = ^ [h^ (R) +/3E, (R)] = j^H^ (R) 



47ry2 
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(12) 



Here y ^ = 1 - j6^ is the relativistic factor of the expanding 
boundary and /oy^^ = (c/47r)//o||.j.. Note that the moving bound- 



ary modifies the surface current which now contains an extra 
factory"^ [3]. 

Consider now briefly the non-relativistic limit of Eqs. (|9]i and 
( fTOb . This limit can be obtained using at ^ ^ and/? —> the 
asymptotic expression ViOITiP) = TiOITifi) = ySVf^ = 
R~ lp~ which yields H~ - //oy, and 

Hp^Ho^coaJl-^Y //^ = -//o^sin^|l + ^J. (13) 

In the lowest order with respect to the factor /? the compo- 
nents of the electric field are given by E^ - /3H^)n(R/p), Er = 
2/3H^)± sin (p{R/p). It is seen that the parallel component H- 
of the magnetic field remains unchanged in the case of non- 
relativistic expansion. 



3. Analysis of the energy balance 

Previously significant attention has been paid ll8l 4l0lll3ill4ll 
to the question of what fraction of energy is emitted and lost 
in the form of electromagnetic pulse propagating outward of 
the expanding plasma. In this section we consider the energy 
balance during the plasma cylinder expansion in the presence 
of the homogeneous magnetic field. When the plasma cylinder 
of the zero initial radius is created at f = and starts expand- 
ing, external magnetic field Ho is perturbed by the electromag- 
netic pulse, H'(r, t) - H(r, f) -Ho, E(r, t), propagating outward 
with the speed of light. The tail of this pulse coincides with the 
moving plasma boundary p - R{t) while the leading edge is at 
p = ct. Ahead of the leading edge, the magnetic field is not 
perturbed and equals Hq while the electric field is zero. 

Our starting point is the energy balance equation (Poynting 
equation) 

where S = ^ [E x H] is the Poynting vector and j = j^e^ + j^e^ 
(with |e^| - le^-l = 1) is the surface current density. The energy 
emitted to infinity is measured as a Poynting vector integrated 
over time and over the lateral surface 5 c of the cylinder with ra- 
dius pc, length /c and the volume Q^. (control cylinder) enclosing 
the plasma cylinder (p^ > R or Q ^ t < pdv). Integrating over 
time and over the volume Q^ Eq. ( fT4l l can be represented as 



Ws(0 = Wy(0 + AWEM(0, 



(15) 



where 



Ws (t) = IcP, 
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(16) 

Here Sp - -^{E^H, - E,H^) is the radial component of the 
Poynting vector. Note that the total flux of the energy over the 
bases of the control cylinder determined by the Poynting's vec- 
tor component 5; vanishes due to the symmetry reason. Wem(0 
and AWem(0 = VVem(O) - Wem(0 are the total electromagnetic 
energy and its change (with minus sign) in a volume Q.c, re- 
spectively. Wj{t) is the energy transferred from plasma cylinder 



to electromagnetic field and is the mechanical work with mi- 
nus sign performed by the plasma on the external electromag- 
netic pressure. At f = the electromagnetic fields are given 
by H(r, t) - Ho and E(r, t) - 0. Hence Wem(O) is the total 
energy of the magnetic field in a volume Q^ and is given by 
W'em(O) = Q - np^lc(HQl%n). Then the change of the elec- 
tromagnetic energy AWem(0 in a volume Q.c can be evaluated 
as 



AW^EM(0 = e- 
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-dr. 



(17) 



In Eq. ( fTTl i Q.'^ is the volume of the control cylinder excluding 
the volume of the plasma cylinder (we take into account that 
H(r, f) - E(r, f) = in a plasma cylinder). Hence the total 
energy flux Ws if) given by Eq. ( fTSI l is calculated as a sum of 
the energy loss by plasma due to the external electromagnetic 
pressure and the decrease of the electromagnetic energy in a 
control volume Q^. For non-relativistic (y6 <K 1) expansion of a 
one-dimensional plasma slab and for uniform external magnetic 
field Ws - IWj - 2AWem, i-c, approximately the half of the 
outgoing energy is gained from the plasma, while the other half 
is gained from the magnetic energy ||8|]. In the case of non- 
relativistic expansion of highly-conducting spherical plasma 
with radius R in the uniform magnetic field Ho the outgoing 
energy Ws is distributed between W j and AWem according to 
Wj = l.Sgo and AWem = O.Sgo with Ws = IQo. where 
2o - HqR^ 16 is the magnetic energy escaped from the spher- 
ical plasma volume MI 311 . Therefore in this case the released 
electromagnetic energy is mainly gained from the plasma. 

Consider now each energy component Ws(t), W]{t) and 
A1Vem(0 separately. Ws(t) is calculated from Eq. ( fT6] l. In the 
first expression of Eq. ( fTSI ) the f'-integral must be performed at 
— ^ f ^ t it < ^) since at ^ f' < — the electromagnetic 
pulse does not reach to the control surface yet and Spipc) = 0. 
From Eqs. (|9]l, (fTOl l and (fTSI l we obtain 
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where rj - pdct < 1. In non-relativistic limit using the asymp- 
totic expression T(ri)/T(P) = P^/j]^ - {t/t^.)^ (with t,, - pdv) 
a.t/3 —> and 77 — > 0, from Eq. ( fTSl l we obtain 



Ws(t) 






(19) 



Next, we evaluate the energy loss Wj(t) by the plasma which 
is determined by the surface current density, j. This current has 
two azimuthal and axial components and is localized within thin 
cylindrical skin layer, R-5 < p < R+6 with (5 — > 0, near plasma 
boundary. Therefore in Eq. ( fTSl l the volume Q^. can be replaced 
by the volume Q^ ~ l^RS which includes the space between 
the cylinders with p - R - 6 and p - R + 6. The surface 
current density is calculated from the Maxwell's equation and 
has been determined in previous section, see Eqs. ( fTTT l and ( fT2l l. 
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Figure 1: (Color online.) The ratios Vs(f) (solid lines) and Yj{{) (dashed lines) for four values of /? as a function of f (in units of t^ 
expressions )18t and J2U with t? = ;r/4. 



Pr/c) calculated from 



As shown below the r-integration of the term j -E in Eq. (fT6t can 
be alternatively expressed via magnetic field. Within the skin 
layer we take into account that E - -i[v x H] and Hp{R) — 0. 
Then 



Qjit)^- r j ■ Efl'r = -^ r V [H X (V X n)]dr 
1 r dE^ r //2 
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where 5^; is the lateral surface of the expanding cylindrical 
plasma. As mentioned in Sec. |2] the factor y"^ in Eq. ( l20l l ap- 
pears because of the moving boundary. In this expression the 
term with ' g^' has been transformed to the surface integral 
using the fact that the boundary of the volume Qg moves with 
a constant velocity v and the electrical field has a jump across 
the plasma surface. Equation ( l20l l shows that the energy loss 
by the plasma per unit time is equal to the work performed by 
the plasma on the external electromagnetic pressure. This ex- 
ternal pressure is formed by the difference between magnetic 
and electric pressures, i.e., the induced electric field tends to 
decrease the force acting on the expanding plasma surface. The 
total energy loss by the plasma cylinder is calculated as 

Wj(t) 1 r^ ,,„ „, f, . Hi] {l-/3^f 
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In a non-relativistic case Eq. ( |2TI ) yields 
Wjit) I 
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(22) 



The change of the electromagnetic energy in a control cylin- 
der is calculated from Eq. ( [TtI i. Ai R < pc < ct (the electromag- 
netic pulse fills the whole control cylinder) we obtain 
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Comparing Eqs. ( fTsT l. dlTT) and ( l23b we conclude that A1Vem(0+ 
Wj{t) = Ws (t) as predicted by the energy balance equation ( fTST l. 
The non-relativistic limit of Eq. (|23] ) can be evaluated from 
Eqs. (dill and ^ using the relation AWEM(f) = ^^^(0 - Wj(t). 
As an example in Fig. [1] we show the results of model calcula- 
tions for the ratios Es (f) = Ws(t)/Qo(t) andr j(t) = Wj{t)/Qo(t) 
as a function of time {pc/c ^ t ^ Pc/v). Here Qoit) = 
nlcR^Hy^n is the magnetic energy escaped from the plasma 
cylinder at time t. For the relativistic factor/? we have chosen a 
wide range of values. We recall that at $J f ^ pdc, i-e. the 
electromagnetic pulse does not yet reach to the surface of the 
control cylinder, Ws (f) - 0. Unlike the case with uniform mag- 
netic field discussed above (see also JSl llSIl ) there are no sim- 
ple relations between the energy components Ws (t), Wj(t) and 



Qo{t). From Eq. (1211 1 it is seen that the ratio rj{t) is constant and 
is given by Tj(t) = (1 + hIJHI)C, where C = [y^li^T{l3)Y^ 
is some kinematic factor For non-relativistic expansion with 
j0 « 1, this factor is C - 1 while at /? ~ 1 this factor is very 
large and behaves as C - (9/8)(l -y6)"' » 1. As expected 
the total energy flux, Ws if), increases monotonically with t. At 
the final stage (f = pdv) of relativistic expansion (withy? ~ 1) 
VKs - W'y. Hence in this case the emitted energy Ws is mainly 
gained from the plasma cyUnder. 
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4. Conclusion 

An exact solution of the uniform radial expansion of a neu- 
tral, infinitely conducting plasma cylinder in the presence of 
a homogeneous magnetic field has been obtained. The elec- 
tromagnetic fields are derived by using the appropriate bound- 
ary and initial conditions, Eq. (|6]l. It has been shown that the 
electromagnetic fields are perturbed only within the domain 
extending from the surface of the expanding plasma cylinder 
p - R - vt Xo the surface of the expanding information cylin- 
der p - ct. External to the cylinder p - ct the magnetic field 
is not perturbed. In the course of this study we have also con- 
sidered the energy balance during the plasma cylinder expan- 
sion. The model calculations show that the emitted energy is 
mainly gained from the plasma cylinder. For relativistic expan- 
sion Ws - Wj and the emitted energy is practically gained only 
from plasma cylinder. 

We expect our theoretical findings to be useful in experi- 
mental investigations as well as in numerical simulations of the 
plasma expansion into ambient magnetic field. One of the im- 
provements of our model will be to include the eff'ect of the 
deceleration of the plasma cylinder as well as the derivation of 
the dynamical equation for the surface deformation. A study of 
this and other aspects will be reported elsewhere. 
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